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Abstract 

The geometry of four-qubit entanglement is investigated. We replace some of the polynomial 
invariants for four-qubits introduced recently by new ones of direct geometrical meaning. It is shown 
that these invariants describe four points, six lines and four planes in complex projective space CP 3 . 
For the generic entanglement class of stochastic local operations and classical communication they 
take a very simple form related to the elementary symmetric polynomials in four complex variables. 
Moreover, their magnitudes are entanglement monotones that fit nicely into the geometric set of 
n-qubit ones related to Grassmannians of Z-planes found recently. We also show that in terms 
of these invariants the hyperdeterminant of order 24 in the four-qubit amplitudes takes a more 
instructive form than the previously published expressions available in the literature. Finally in 
order to understand two, three and four-qubit entanglement in geometric terms we propose a 
unified setting based on CP 3 furnished with a fixed quadric. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.65.Ta, 02.40.-k 
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I. INTRODUCTION 



Recently the problem of characterizing n-qubit entanglement classes has generated con- 
siderable interest. This problem was raised within the context of quantum information 
theory regarding the physical phenomenon of entanglement as a resource. In order to ex- 
ploit the capabilities encoded in this resource for different tasks of quantum information 
processing we have to somehow measure it. During the past few years a number of useful 
entanglement measures for pure states has been appeare d 1 ! 2 ! 3 ^ 5 ! 6 ^ 8 ' 9 . Classifications up to 
four qubits has been appeare d 10 i 11 i 12 and the interesting geometric structures associated with 
entangled pure states has been note d 9 i 12 i 13 i 14 i 15 i 16 i 17 i 18 . 

In our previous set of pape ra 9 ' 16 ' 17 using some results from twistor theory we initiated 
an approach for understanding n-qubit entanglement in geometric terms. We have shown 
that this problem can be completely solved for three-qubits, and we obtained partial results 
for a special subclass characterized by n-qubit entanglement monotones. The aim of the 
present paper is to add some interesting new results on the geometry of four-qubit SLOCC 
(stochastic local operations and classical communication^) invariants. Such invariants has 
been introduced in [5], where the Hilbert series of the algebra of invariants has been found. 
This result enabled the authors to construct a complete set of four algebraically indepen- 
dent invariants of degree 2,4,4,6 in the complex coefficients characterizing the four-qubit 
entangled pure state. The values of these invariants on the SLOCC orbits of Ref. [11] were 
also given. Moreover, an explicit formula for the hyperdeterminant of degree 24 was also 
obtained. The authors of Ref. [5] has conjectured that some of the invariants might have a 
geometric meaning. In this paper we show that this is indeed the case. 

In Section II. we present a new set of four invariants, by replacing two from the ones 
of Ref. [5]. In Section III. we clarify the geometric meaning of these invariants in terms 
of the geometry of CP 3 the complex projective space. In Section IV. we show that using 
our new set of invariants the expression for the 24th order hyperdeterminant takes a more 
instructive form than the one that can be found in [5]. Moreover, it turns out that the 
entanglement classes invariant under SLOCC transformations take a very simple form related 
to the elementary symmetric polynomials in four complex variables. The magnitudes of our 
invariants turn out to be entanglement monotones that fit nicely into the geometric set of 
n-qubit entanglement monotones related to Grassmannians of /-planes in C L with L = 2 n , 
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I < L found recently. Finally our conclusions and some comments are left for Section V. 



II. INVARIANTS 



Let us write an arbitrary four qubit state in the form 



|*>= Z l]kl \i 3 kl) e C 2 ® C 2 ® C 2 ® C 2 (1) 

i,j,k,l=0 

where \ijkl) = <S> \j) 2 ® 1^)3 ® 104- Following [5] we introduce decimal notation for 
Zijki = Z r where r = 8i + 4j + 2k + I and the matrices 



C 



= (a, b, c, d) 



(2) 



M 




(3) 



^ Zq Z4 Z$ Zi 2 \ 

Z\ Z 5 Zc, Zi 3 

Z 2 Zq Z 10 Z14 

\Z 3 Z 7 Z n Z 15 J 

( Zq Z 2 Z$ Z\q\ 

Z\ Z 3 Z§ Z\\ 

Z^ Z§ Z\ 2 Z14 

\Z 5 Z 7 Z 13 Zi 5 J 

^ Zq Z\ Z$ Zg^ 

Z 2 Z 3 Z w Zu 

z 4 z 5 z 12 z 13 

\Z% Z7 Z\4 Z15J 

Here A, B, C, D e C 4 are considered as four column vectors and A,B,C,V are 2x2 
matrices with T referring to transposition. We wish to describe the geometry of four qubit 
entanglement in terms of the four vectors A a , Bp, C 1 and D$ living in C 4 where a, (3, 7, 5 — 
0, 1,2,3. Hence we regard the matrix C as fundamental. The matrices M and Af will be 
used later. 

We are interested in studying a subset of polynomials in the complex numbers Z r ,r = 
0, ... 15 that are invariant under the SLOCC group of stochastic local operations and classical 




(4) 
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communication i.e. SL(2, C)® 4 . Such transformations are of the form 

\*)*-+{S 1 ®S 2 ®S 3 ®S i )\V), (5) 

where S m G SL(2,C),m = 1,2,3,4, and takes the (JTJ) form with the indices of S m 
referring to the label of the C 2 in the tensor product they are acting on. 

In order to define the SLOCC invariants we wish to propose we will introduce two extra 
structures on C 4 . The first one is a bilinear form g : C 4 x C 4 — > C such that for two vectors 
A, B G C 4 we have 



(A, B) ^ g{A, B) = A B = g aP A a B p = A a B a 



(6) 
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to A 




U A 




M 


l-l oJ 


U o; 







(7) 



where 

/o l\ 
0-10 
0-100 
\1 0/ 

a, (3 = 0, 1,2,3 and summation for repeated indices is understood. Notice that the matrix 
of our symmetric bilinear form can be written in a tensor product form g = e <8> s where e is 
invariant under the group SX(2,C), i.e. we have SeS T = e with S G SL(2, C). Moreover, 
since 



A a — Z 00 ki, Bp — Z m ki, C 1 = Zioki, Ds — Znki, a, /3, 7, S = 0, 1, 2, 3, k,l — 0,l 

(8) 

quantities involving this symmetric bilinear form are automatically invariant with respect 
to SL(2, C) <S> SL(2, C) transformations of the the third and fourth qubit, i.e. the ones of 
the form / ® I ® S 3 <8> S 4 . Since SL(2, C) <8> SL(2, C))/Z 2 ~ SO(4, C) it follows that greek 
indices like a, (3 etc can also be regarded as vector indices under 5*0(4, C). This conversion of 
complex four- vectors into complex 2x2 matrices has already been used elsewhere to connect 
the results of twistor theory to the geometry of entanglement*^. Hence the columns of the 
matrix L of (0) transform as vectors under transformations of the form I ® I ® S3 ® £4 and 
as the 00, 01, 10 and 11 components of a tensor under the ones with the form Si® £2® I® I- 
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The second structure, as we will see, is related to the notion of duality in CP 3 . For the 
vectors A, B, C and D let us introduce their duals as 

a a = —ea^sB^C^D 5 , bp = e a i3 1 sA a C 1 D 6 , c 7 = e a ^A a B 13 ' D & , d$ = — e aj 3 1 sA a B l3 C 1 . 

(9) 

Clearly these quantities transform as vectors under transformations of the form 7® 7® 63(8) £4 
and a straightforward calculation shows that a, b, c and d behave under the ones of the form 
Si ® S 2 ® 7 ® 7 as the 11, 10, 01 and 00 components of a tensor respectively. 

Let us now introduce the notation A AB (a bivector) corresponding to the antisymmetric 
matrix A a Bp — ApB a and the one 

L = Det£ = e^AaBpC^Ds. (10) 
Now the SLOCC invariants we wish to propose are 

/i = i(A-D-B-C). (11) 

h = l - ((A A B) • (C A D) + (A A C) • (B A D) - 1(A A D) 2 -^(BA C) 2 )^) (12) 

7 3 = i(a-d-b.c). (13) 

h = L (14) 
Here quantities like (A A B) • (C A D) are defined as 

(AAB)-(CAD) = (A a Bp-A p B a )(C a D p -C p D a ) = 2((A-C)(B-D)-(A-D)(B-C)). (15) 
The first of our invariants Ii takes the form 

h = -H = -(Z Zi 5 — ZiZ u — Z 2 Zi 3 + Z 3 Z U — Z4Z11 + Z 5 Z W + Z e Z 9 — Z 7 Z 8 ), (16) 
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showing that 1\ = \H where H is one of the basic invariants of Ref. [5]. It is just a special 
case of the n-tangle with n even introduced earlier by Wong and Christensen^. Reverting 
to binary notation it is easy to show that I\ is also a permutation invariant^**^. 

Our last invariant 1^ is just L = Det£ an invariant also introduced by the authors of [5]. 
These authors have also introduced two more invariants of order four and six denoted by M 
and D respectively (the first of them being just minus the determinant of our matrix Ai. 
They have shown after obtaining the Hilbert series that the invariants H, L, M and D are 
algebraically independent and complete. Here instead of the invariants M and D we prefer 
the new ones I2 and I3. As we will see the set (ii, I2, I3, 1 a) is of geometrical significance. 
Moreover, it turns out that the values of these invariants on the generic SLOCC orbit— of 
four-qubit entangled states are just the elementary symmetric polynomials in four complex 
variables. We will also show that in terms of this new set of invariants the explicit formula 
for the hyperdeterminant of degree 24 takes a more instructive form then the corresponding 
one of Ref. [5] . 

III. THE GEOMETRIC MEANING OF FOUR QUBIT INVARIANTS 
A. The invariant I2 

Let us now explain the structure of ^ In order to do this we introduce another C 4 
corresponding to the four- vector structure also present in the first two indices of Z^. Con- 
verting the first two spinor indices to vector ones labelled by \i = 0, 1, 2, 3, what we obtain is 
a "vector- valued" four vector Z m = (A, B, C, D) T . (Alternatively regarding Z m as a 4 x 4 
matrix we obtain the matrix C of Eq. (2).) Let us now also supply this new copy of C 4 
with the bilinear form g known from Eq. (jHJ with matrix /i, v — 0, 1, 2, 3. Define now 
the second exterior power of a matrix as the map 

yy 2 . C nxn _^ c (2)x(£) ( 17 ) 
which takes an M^ u e C nxn , < fM, v <n— lto 

= {f\M) Lj = M^M^ - M^M^ (18) 
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where I = {//i, 1x2} with 0</ii</i2<^ — 1 and J = {ui, z/ 2 } with < v\ < v% < n — 1. 
For the 4x4 matrix of our bilinear form g we have 



(2)IJ 



(0 














A 














1 











-1 




















-1 











1 














v 1 














0/ 



(19) 



where J, J = 01, 02, 03, 12, 13, 23 

Let us now introduce the " bivector- valued" 4x4 antisymmetrix Pliicker matrix 

/ A A B AACAAD\ 
-A A B BAC BAD 

-A AC -B AC CAD 

\-AAD -BAD -CAD j 



P 



(20) 



Notice that since the elements of P„ v are separable bivectors it has the index structure P, 



Explicitly we have 



fiua/3 — Z m Z v fi — Z^Z va . (21) 

They are the Pliicker coordinates of six lines in CP 3 labelled by the subscripts jiu = 
01,02,03,12,13,23. Alternatively one can introduce the " bivector- valued" vector P with 
six components Pi, I = 01, 02, 03, 12, 13, 23 with the same geometric meaning. 



Pi = (A A B, A A C, A A D, B A C, B A D, C A D) T . 
In this notation the invariant I2 of Eq. (12) can be written in the form 



(22) 



G IJ Pi ■ Pj 



6 ^ 



_pixvo.fi p 



[j,va(3- 



(23) 



12 u 6 6 

Notice that the elements of the Pliicker matrix being now separable bivectors, are multiplied 
together according to the (fTK|) law bringing in the • product. Since I2 contains contractions 
with respect to four SL(2, C) ® SL(2, C) invariant matrices g (two of them operates on 
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the first two and the other two on the last two qubits) it is automatically an SL(2, C)® 4 
invariant. 

From the form of the Pliicker matrix Eq. (20) it is also clear that it encapsulates infor- 
mation concerning four qubit entanglement in terms of six planes in C 4 or alternatively six 
lines in CP 3 . Hence we managed to identify the SLOCC invariant I 2 as a line invariant in 
complex projective space. 

Let us now clarify the relationship of our invariant I 2 with the ones of Ref . [5] . There 
also the invariants L = Det£, M = DetM and N = DetA/" were defined where the 4x4 
matrices are given by Eqs. (2-4). Notice that our convention for M differs in a sign from 
the one adopted in Ref. [5]. It can be shown that M = L + N as can be verified by 
calculating the determinants of the matrices in Eqs (3) and (4) containing 2x2 blocks. Now 
a straightforward calculation shows that in terms of the algebraically independent invariants 
H — 27i, L — I4 and M preferred by ref. [5] we have the relation 

6I 2 = H 2 + 2L - AM. (24) 

B. The invariant ^3 

In order to understand the geometric meaning of the invariant I3 of Eq. (13) we have to 
recall some results from projective geometry. A plane in CP 3 consists of a set of points with 
homogeneous coordinates X a a — 0, 1, 2, 3 which satisfy a single linear equation of the form 

a a X a = (25) 

where the complex numbers a a are called the coordinates of the plane. Clearly a a and \a a 
with A 7^ determine the same plane in CP 3 so the set of planes in CP 3 is itself a CP 3 called 
the dual projective space. A plane in CP 3 is a CP 2 . There is a unique plane containing 3 
general points in CP 3 . If Bp, C 7 and D$ are 3 general points then there is a unique solution 
up to proportionality of the three equations 

a a B a = a a C a = a a D a = (26) 

given by the first of Eq. (9). Now we see that the four- vectors a a , bp, c 7 and d$ defined in 
Eq. (9) are the (dual Pliicker) coordinates of four planes in CP 3 . They are defined by the 
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three points (B,C,D), (A, C,D), (A, B,D) and (A, B, C) respectively. Alternatively for 
these planes we can use the Pliicker coordinates (separable trivectors) BACAD, AaCAD, 
AABAD, AABAC where for example 

(AABAC) a/ 3 7 = AaBpCy+AjBaCp+ApB^Ca— ApB a C 7 — A^BpCu— A a B 7 Cp = 3\A[ a BpC 7 ]. 

(27) 

The set of points common to two planes is a line in CP 3 . A line is given by the points 
satisfying two linear equations of the form 

a a X a = b a X a = 0. (28) 
A sufficient and necessary condition for these equations to hold is 

a [a b^ = 0. (29) 

Hence in order to characterize the line (a CP 1 ) defined by this equation we can either use 
the dual Pliicker coordinates (a separable bivector) 

(a A b) a(3 = a a bf3 - a p b a = 2!a [a 6 /3] (30) 
or using Eq.(9) the Pliicker coordinates (another separable bivector) 

(CA-D) l5 = C 1 D s -CsD 1 = 2\C b D 5] . (31) 

This example shows that e.g. the planes with the Pliicker coordinates B A C A D and 
A A C A D intersect in the projective line given by the Pliicker coordinates CAD. 

Now clearly /3 of Eq. (13) is an SL(2, C)® 4 invariant. Indeed, the dual Pliicker coor- 
dinates a, b, c and d are transforming as the 11, 01, 10 and 00 components of a tensor 
under transformations of the form Si <g) S% <S> I <8> I, hence the combination a • d — b ■ c is 
an invariant with respect to such transformations due to the "determinant-like" structure. 
Moreover this quantity is also invariant under transformations of the form I <g) I <g) S3 <S> S4 
due to the occurrence of the SL(2, C)® 2 invariant • product of Eq.(6). Notice that ^3 has 
the same structure as I\. This exemplifies a general pattern: suppose we have an invariant 
(Ji) , then find a set of covariants (i.e. a, b, c and d) to construct a new invariant (I3) by 
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exploiting the existing structure of the original invariant (ii). It is clear that I± describes 
constellations of points and I 3 describes planes in CP 3 . It is important to realize however, 
that unlike ^3 the invariant I\ is also a permutation invariant. 

Let us also express our invariant in terms of the Pliicker coordinates A A B A C etc as 

I 3 = ((A A C A D) • (A A B A D) - (B A C A D) • (A A B A C)) . (32) 

Alternatively like in Eq. (20) one can define a third order totally antisymmetric Pliicker 
tensor P^ up with trivectors as elements. The four independent elements of P^p are P i2 = 
A A B A C, P123 = B A C A D, P 023 = A A C A D and P 013 = A A B A D. Notice that these 
quantities have the index structure -£0120/37 = (A A B A C) Q| g 7 where the definition of Eq. 
(27) holds. Using this notation we have 

1 1 

~ 12 2 VP ~ 12 2 r V.vpaP~fi \°°) 

an expression to be compared with Eq. (23) obtained for our line invariant. 

Finally let us relate our invariant ^3 to the sixth order ones of Ref. [5]. Define the 
quadrilinear form 

1 

Z(x,y,z,t)= ^ ZijkiXiVjZkti- (34) 

i,j,k,l=0 

Using this for each pair of variables one then defines the covariants (now we define it for the 
pair xy) 



b xy (ic,z) = Det ^ 



d 2 Z 

dzidtj 



(35) 



Now let us reinterpret these biquadratic forms as bilinear forms on S 2 C (the symmetric part 
of C 2 ® C 2 ) as 



yoyi 



(36) 



b xy (x,y) = (x ,x x 1 ,x 1 )B xy 
i.e. B xy is a 3 x 3 matrix. Then following Ref. [5] we define 

D uv = Det(B uv ). (37) 
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Hence we have six sextic invariants D xy , D zt , D xz , D yt , D xt and D yz . According to Ref. [5] 
only four of them is independent due to the relations D xy = D zt , D xz = D yt and D xt = D yz . 
Now a straightforward calculation shows that 

h = ^{D xz + D xt ). (38) 

In Ref. [5] the authors used the invariant D = D xt as a fundamental one satisfying the 
relation D xz — D xt = HL hence we can write 

h = D+ l -EL. (39) 

By virtue of Eqs. (14), (16), (24) and (39) the relationship between our set of invariants 
(ii, I 2 , h, h) and the ones (H, M, L, D) used in Ref. [5] is established. 



C. The invariants I4 and I2 

The meaning of this invariant is clear. L = DetC is vanishing when the vectors A, B, C 
and D are linearly dependent. Since I4 and I2 are both of fourth order let us now explore 
the relationship between them. 

Let us label the six lines as in Eq. (20). Hence for example P i is the line corresponding 
to the separable bivector AAB. More precisely this object has the index structure (Poi) a p — 
(A A B) a/3 = A a Bp — ApB a . In this notation 

L = \e a ^ s (P 01 U(P 23 ), 5 . (40) 

We can regard this expression as a symmetric bilinear form in the six Pliicker coordinates 
of the two lines A A B and C A D. Let us denote this bilinear form by (, ) hence we have 

[P fJ,v 1 P per) 1 ^ ( Pfiv j Ppcr ) • 

(41) 

Defining the dual of a bivector as 

*P a( s = \e a p lS p-< 5 , (42) 

it is easy to show that 

{P^,P p(T ) = {*P^*P pa ). (43) 
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In this notation the equation (P pu , P up ) = expresses the fact that the planes described by 
the separable bivectors P^ v and P up in C 4 have a line in common. In the CP 3 picture this 
is equivalent to the fact that the corresponding lines in CP 3 intersect in a point. Hence we 
can look at L as a line invariant too, moreover we have the obvious relations 

h = (Poi, P23) = (P 2,P 3 l) = (P 3,Pl2> = CP01/P23) = (*P02,*P 3 l) = (*P03,*Pl2). (44) 

Let us now look at quantities like {P pu , *PpaY- We have for example 

(Poi, ^23) = (A • C)(B • D) - (A • D)(B • C = ±(A A B) • (C A D). (45) 
Since (P^, *P pa ) = (*P^u, Ppa) one can write now the invariant QI 2 in the form 

6/2 = (^01,^23) + CPoi,^)) + (P 2,*Pl 3 ) + (*P02,Pl 3 ) ~ (P 03 ,*P 03 ) ~ (Pl2,*Pl2)- (46) 

Using Eq. (43) let us calculate 6(h ± h)\ We get 

6(I 2 ±h) = ±(P 01 ±*P 01 ,P 23 ±*P 23 )±(P 02 T*P 2,P 3 iT*P 3 i)-(Pi2T*P 03 ,*Pi2TP 03 )- (47) 

As we will see these invariants will occur in the expression for the hyperdeterminant. The 
geometric meaning of these invariants is connected to the intersection properties of the self- 
dual (*P = P) or anti- self- dual (*P = —P) parts of the planes in C 4 (or alternatively of 
lines in CP 3 ). For example for Pqi self-dual, P31 anti-self-dual and P12 identical to the dual 
line of P30 (an equivalent condition for this is L(B A C) = — b A c) the invariant I4 — I 2 
vanishes. It is easy to check that the invariant U occurring in Ref. [5] can be related to one 
of these invariants as 

U = H 2 -4(L + M) = 6(/ 4 - h). (48) 

The fact that (among others) this invariant might have a geometric meaning was raised in 
Ref. [5]. 

IV. THE HYPERDETERMINANT 

Let us now consider the hyperdeterminant P 4 for the four-qubit system. As it is well- 
known for two-qubit systems the determinant D 2 = Z 00 Z n — Z iZ w is related to the 
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concurrence 2 as C = 2 1 -Z3>2 1 characterizing two-qubit entanglement. Similarly for three-qubits 
the basic quantity is the three-tangle 2 r = 4|Z^3 1 which is related to the hyperdeterminant 
D 3 of a 2 x 2 x 2 tensor formed from the 8 complex amplitudes Z^. D 3 is an irreducible 
polynomial in the 8 amplitudes which is the sum of 12 terms of degree four. For the explicit 
expression of D 3 see e.g. the book of Gelfand et.al.— . It is known that the next item in 
the line namely the hyperdeterminant of format 2x2x2x2 is a polynomial of degree 
24 in the 16 amplitudes Zyjy which has 2894276 terms^i. An expression in terms of the 
fundamental invariants (H,L,M,D) was given in Ref. [5]. Here we are interested in the 
explicit form of -D4 the hyperdeterminant of the 2x2x2x2 tensor Zijki based on the special 
invariants (ii, I2, 13, I4) we have found in our CP 3 picture. 

As it is well known^ the hyperdeterminant -D4 is the unique irreducible polynomial in 
the 16 unknowns Zijki that vanishes whenever the system of equations 



F — Z 00 q + Z oom t + Zqq W z + Z i 00 y + Zi 000 x + Z oon zt + Z mm yt + • • • + Z nw xyz + Zmixyzt 



has a solution (x , yo, z , t ) in C 4 . 

Using the method of SchlafTi according to Theorem 14.4.1 and Corollary 14.2.10 of Ref. 
[20] -D4 coincides with the discriminant A of D^iZ^i + ^^ijfc«) considered as a polynomial 
in A divided by 256. This method has already been used to obtain a much simpler form for 
D 3 of geometric meaning^. For _D 4 a method equivalent to this has been applied with the 
result 5 




dF _dF _dF _ OF 

dx dy dz dt 



(49) 



where 



(50) 



256£> 4 = S 3 - 27T 2 



(51) 



where 



125 = U 2 - 2V, 216T = U 3 - 3UV + 216£> 2 , 



(52) 



with 



U = H 2 -4(L + M), V = 12(HD + 2LM). 



(53) 
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Let us now express D4 in terms of our invariants (Ii, J 2 , I3, I4) of geometric significance. 
Using relations (14), (16), (24) and (39) in Eqs. (51-53) we obtain the result 

S = {II - P 2 ) + 4(/ 2 2 - hl 3 ), T = {II - P 2 ){lf - I 2 ) + (I s - hl 2 f. (54) 

In this form it is obvious that the combined invariants I 4 ± / 2 , If — I2 of fourth (I3 — I1I2) 
of sixth and If — I1I3 of eight order should play a basic geometric role. We have already 
clarified the geometric meaning of the first two invariants. They are related to self-duality 
and anti-self-duality of the corresponding lines in CP 3 . One of these invariants I 4 — I2 is 
just |{7 also used in Ref. [5]. In our form of D4 we prefer to also use the dual combination 
I4 + 1 2 . For the time being we do not know any geometrical interpretation of the other 
combinations. Intuitively it is clear that the invariant If — 12 should play a similar role than 
the other fourth order invariants. Indeed we have chosen the third and fourth qubits to play 
a special role. (An equivalent picture arises when a special role is assigned to the first and 
the second qubit.) The C 4 defined by them is equipped with the bilinear form Eq. (6). The 
null vectors (i.e. the ones satisfying X ■ X = 0) describe a quadric embedded in CP 3 which 
is isomorphic to CP 1 x CP 1 , i.e. it is ruled by two families of projective lines which can be 
shown to be self-dual or anti-self dual respectively^. Projective lines lying entirely inside a 
fixed quadric are called isotropic lines. Had we chosen the two qubits playing a special role 
differently the notion of self or anti-self-duality of isotropic lines would have been defined 
with respect to a different quadric. In this picture we conjecture that the invariants I4 — I2 
and If — I 2 would play a dual role. Since altogether we have three inequivalent choices 
then we can conclude that the fourth order invariants are related to the notion of duality 
of isotropic lines with respect to a fixed quadric in CP 3 . It would be interesting to find a 
geometric interpretation for the remaining invariants too. 

Let us now consider another interesting property of D4 expressed in terms of our invariants 
(Ji, hyh, I±)- As it is well-known the discriminant A of the polynomial e^w 4 + e 3 u> 3 + e2W 2 + 
t\W + e is given by the expression (see Ref. [20] Eq. (1.35) on page 405) 

A(e 4 w + e 3 w + e 2 w + e\W + e ) = 256ege 4 — 192e eie 3 e 4 — 128e e 2 e 4 

00 0/1 00 00 o 

+ 144e e 2 e 3 e4 — 27e e 3 + 144e e 1 e 2 e 4 — 6e e 1 e 3 e 4 — 80e eie 2 e 3 e4 
+ 18e eie 2 e 3 + 16e e 2 e 4 - 4e e 2 e 3 - 27e 4 e 4 + 18e 3 e 2 e 3 e 4 - 4e 3 e 3 
- 4e 2 e 2 e 4 + e\e\e 4 . (55) 
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Let us now consider the polynomial of the special form 

p[I x , I 2 , h, h- w}=w A - (Ah)w 3 + (QhW - (4I 3 )w + l\ (56) 
Then a straightforward calculation shows that 

256D 4 = A(p[I 1 ,I 2 ,I 3 ,I 4 ;w)). (57) 

Notice that the polynomial p is not directly related to the one arising from the method 
of Schlafli. In this case one obtains 

D 3 (Z 0jkl + XZ ljkl ) = h 4 X 4 + M 3 + M 2 + hX + h (58) 

where unlike the ones e s the coefficients h s s = 0, 1 ... 4 are fourth order polynomials 
of the Zijki that are not invariant with respect to the full group SL(2, C)® 4 . However, the 
discriminant of this polynomial again gives 256-D4 which is already an invariant with respect 
to the full group of SLOCC transformations. 

In order to illustrate the advantages of using our invariants (Ji, J 2 , J 3 , J 4 ) let us now 
calculate their values on the generic SLOCC class. A generic pure state of four qubits can 
always be transformed to the formi 1 - 

\G abcd ) = i±^(|0O0O> + |llll» + ^(|0Oll> + |110O» 

+ ^(|0101) + 11010)) + ^(|0110) + 10110)), (59) 

where a, b, c, d are complex numbers. For this state the reduced density matrices obtained 
by tracing out all but one party are proportional to the identity. This is the state with 
maximal four-partite entanglement. Another interesting property of this state is that it 
does not contain true three-partite entanglement 1 *" A straightforward calculation shows 
that the values of our invariants (Ji, I2, h, h) occurring for the state \G a bcd) representing 
the generic SLOCC class are 



I, = \y + b 2 + c 2 + d 2 ], h = \[{abf + {acf + (ad) 2 + (be) 2 + (bd) 2 + (cd) 2 ], (60) 
4 6 



J 3 = -[(abc) 2 + (abd) 2 + (acd) 2 + (bed) 2 }, h = abed, (61) 
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hence the values of the invariants (4Ji, 6J2, 4/3, J|) occurring in the polynomial Eq.(56) are 
given in terms of the elementary symmetric polynomials in the variables (xi, x%, xs, x±) = 
(a 2 , b 2 , c 2 , d 2 ). From this and Eq. (57) it immediately follows that the value of the hyperde- 
terminant on the SLOCC orbit represented by the state |G a 6cd) is 

D * = 2o6 ni<j(Xl " Xj)2 = 2o6 V(a2 > ^ ° 2 ' ^ )2 ' (62) 
in accordance with Ref. [5] where V is the Vandermonde determinant. For the other SLOCC 
classes and the values of the invariants (H,L,M,D) see Ref. [5]. It is straightforward to 
give the alternative values of (Ji, J 2 , I3, 1 4) on these classes. 

V. COMMENTS AND CONCLUSIONS 

In this paper we have considered some aspects of the problem of understanding four 
qubit entanglement in geometric terms. We have replaced two from the set containing four 
algebraically independent invariants of Ref. [5] by new ones. In this way all four invariants 
have a simple geometric meaning. I\ is based on 4 0-planes (points), I2 on 6 1-planes (lines), 
I3 on 4 2-planes (planes) and finally I4 on a single 3-plane in CP 3 . According to Theorem 
2. of Ref. [24] the magnitudes of these invariants can be used as entanglement monotones 
characterizing four-qubit entanglement. Moreover, for an arbitrary four-qubit state after 
calculating the set of invariants (ii, I2, 13, I4) and the value of the hyperdeterminant D4 , 
for D 4 7^ we obtain four different roots of the fourth order equation Eq. (56). These 
roots are just the complex numbers (a 2 , b 2 , c 2 , d 2 ) . Their square roots produce the values 
(±a, ±6, ±c, ±0?) appearing in the canonical form G a b c d- This shows that the study of the 
degenerate cases of multiple roots of Eq. (56) arising for _D 4 = 0, could be useful for obtaining 
the parameters of the canonical formal from the values of the basic invariants. This process 
is similar to the spirit of the one found for the three-qubit caseS. There by calculating the 
values of the independent SU(2)® 3 invariants the canonical form of an arbitrary three-qubit 
state was found. 

We note that our four qubit entanglement monotones fit nicely into the scheme of Ref. 
[9] generating a class of n-qubit entangled monotones based on bipartite decompositions of 
Ji = C 2 ". The basic idea followed there was to consider the manifold of subspaces of H i.e. 
suitable Grassmannians with the corresponding Plucker coordinates for them. Equivalently 
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we should consider subspaces of the corresponding projective spaces P(TC). Fixing a quadric 
Q defined by a bilinear form similar to Eq. (6) a class of SLOCC invariants expressed in 
terms of these Pliicker coordinates can be generated. In this way we were able to reproduce 
three of the basic four-qubit invariants (i.e. the triple (H, L, M)). Now we see that by 
employing also the notion of projective duality all four algebraically independent invariants 
of the four-qubit case can be written in the Pliicker form. (See Eqs. (33) and (23) for our 
new invariants). 

Finally let us propose a suggestive geometric picture for four-qubit entanglement. In the 
usual picture^ 2 - a four-qubit state can be represented by a single point in CP 15 . Different 
SLOCC classes correspond to this point lying on different subvarieties in CP 15 . Here we 
would like to suggest an alternative picture. To a four-qubit state we associate a set of four 
points A, B, C, D, six lines AAB, AAC, AAD, BAC, BAD, CAD and four planes 
AABAC, AACAD, AABAD, BACADinthe space CP 3 of smaller dimension. It is 
easy to see looking at the intersection properties of these geometrical objects that they are 
forming a tetrahedron in CP 3 . This correspondence between entangled states and geometric 
objects (unlike the previous one) is nonlocal. The invariants (ii, I2, 13, 1 a) we have proposed 
obviously characterize the properties of this terahedron. For example for J 4 = the four 
points corresponding to the four vectors A, B, C and D in C 4 are not linearly independent 
(some of them are proportional), hence the tetrahedron is degenerating to a triangle or a 
line etc. depending on the degree of degeneracy. We conjecture that the vanishing of the 
other three invariants somehow characterize more intricate degeneracies occurring with lines 
and planes of the tetrahedron. The class containing no degeneracy is the G^ed-class that 
can be represented by a regular terahedron. It would be interesting to understand how the 
SLOCC classes arise in this picture. 

Let us elucidate the meaning of the proposed correspondence a little bit further. The 
usual geometric classification schemes for multiqubit systems are based on the use of hy- 
perdeterminants of more general type. These hyperdeterminants describe geometrically 
hypersurfaces projectively dual to the so called Segre embeddin g 12 ' 21 ^ representing the sub- 
variety of totally separable states 1 ^. For two, three and four-qubits for instance the manifold 
of totally separable states is CP 1 x CP 1 embedded in CP 3 jCP 1 x CP 1 x CP 1 in CP 7 , 
and CP 1 x CP 1 x CP 1 x CP 1 in CP 15 respectively. Here n = 2,3,4 qubit states carrying 
entanglement are represented by points off the Segre surfaces. Such surfaces are representing 
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totally separable states in different projective spaces (CP 2 Here following the spirit of 
our previous set of pape ra 9 i 16 i 17 i 22 we prefer to suggest a unified CP 3 picture. For n = 2, 3, 4 
we take CP 3 with a fixed quadric Q based on our choice of bilinear form Eq. (6). Points 
lying on Q as a subvariety of CP 3 correspond to null vectors in C 4 . For n = 2 to a separable 
or an entangled state corresponds a point on or off Q respectively. For n = 3 we get the 
following geometric picture^^. To a three-qubit state in the GH Z-cl&ss corresponds a line 
in CP 3 intersecting Q at two points. To a state in the iy-class corresponds a line tangent 
to Q at a point. The separable classes B(AC) and C(AB) are represented by isotropic lines 
lying entirely in Q. They are self-dual and anti-self-dual lines belonging to the two different 
rulings of Q. The A(BC) and (A)(B)(C) classes again correspond to the degenerate case of 
points on and oif the quadric Q (see also the pictorial representation of Ref. [22]). We expect 
a similar pattern to exist also for the four-qubit (n = 4) case. Here we have more lines 
arranged to form a tetrahedron, and we have to consider constellations of these lines with 
respect to our fixed quadric Q. The picture arising in this way has some striking similarity 
with the Majorana 23 representation of states with spin s. One can represent geometrically 
a state of spin s as a single point in CP 2s , or alternatively as a constellation of 2s points on 
CP 1 i.e. the Bloch sphere. Some degeneracies can occur in this case e.g. when 2s points 
degenerate to a single one with multiplicity 2s corresponding to the states of highest and 
lowest weights. In the same spirit we would rather represent n-qubit entangled states in 
CP 3 . Here, however in order to account for the nonlocality of multiqubit quantum entan- 
glement instead of merely a collection points we have to also consider constellations of lines 
and planes in CP 3 . Though this analogy is very appealing we expect it to run out of steam 
for the n = 5 case (five-qubits) where probably we should furnish CP 3 with more extra 
structures than a quadric. However for n > 5 in principle we can consider constellations of 
simplexes in CP 3 related to fundamental invariants of the SLOCC group whose combina- 
torial variability should somehow correspond to the proliferating number of entanglement 
classes. 
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